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The next-to-next-to-leading order (NNLO) QCD analysis of the revised experimental data of the CCFR col- 
laboration for the xFs structure function of the deep-inelastic scattering of neutrinos and antinuetrinos on the 
nucleons is made by means of the Jacobi polynomial expansion technique. The NNLO values of the QCD coupling 
constant are determined both without and with twist-4 contributions taken into account. Theoretical ambiguities 
of as are fixed using the methods of the Fade approximants. We observe that the [0/2] Fade approximant is more 
appropriate for this purpose than the [1/1] one. 



1. In the last several years the considerable 
progress was achieved in the area of the ana- 
lytical evaluation of the characteristics of DIS 
at the NNLO of perturbative QCD in the MS- 
scheme. Indeed, the NNLO corrections were 
found in Ref. |l| for the coefhcient functions of the 
structure function (SF) F2, while in the works 
of Ref.[^ the similar NNLO calculations were 
made for the SF XF3. The results of ReL|] 
are in agreement with the expressions for the 
NNLO corrections to the coefficient functions of 
the non-singlet (NS) moments of F2 SF, calcu- 
lated analytically in the case of the even moments 
with n = 2,4,6,8,10 ^J. Moreover, in the case 
of n = 2,4,6,8,10 the analytical expression of 
the NNLO corrections to the anomalous dimen- 
sions of the NS moments of F2 are also known 
at present [Q. Together with the expression for 
the NNLO correction to the QCD /?- function [D , 
the resuhs of Ref. §-Ref. g are forming the theo- 
retical background for the extraction of the value 
of the parameter from the NS fits of the 

DIS data at the level of new theoretical precision, 
namely with taking into account the effects of the 



NNLO corrections. 

The classical approach of the incorporation of 
the information about the QCD approximations 
of the SFs into the analysis of the concrete DIS 
data is the DGLAP approach It is known 
that in its standard variant the DGLAP machin- 
ery is based on the solution of the corresponding 
integro-differential equations. However, at the 
present level of the development of the calculation 
technique it is still impossible to use this vari- 
ant of the DGLAP approach at the NNLO level. 
Indeed, in spite of the fact that the NNLO ap- 
proximations of the DIS cocfhcicnt functions are 
already known from the results of Refs.|0,|^, the 
explicit NNLO expressions of the corresponding 
kernels are not yet calculated. This still missed 
information can be obtained only after the ana- 
lytic calculations of the NNLO corrections to the 
anomalous dimensions of the relevant composite 
operators in the case of arbitrary number of the 
moments n will become available. 

However, there are the methods, which are al- 
lowing to overcome this limitation and to take 
into account the results of the NNLO calcula- 
tions of Rcfs.@]-||| in the process of the NNLO 
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analysis of the DIS data within the framework of 
the MeUin image of the DGLAP equation. These 
methods are based on the apphcation of the infor- 
mation about the NNLO renormahzation group 
evolution of the finite number of the Mellin mo- 
ments of the DIS SFs and on the further approx- 
imate reconstruction of the behavior of the SFs 
themselves. The first realization of the idea of 
the reconstruction of the SFs from the sets of 
the finite number of their Mellin moments was 
proposed in Ref. 0| , where the series in Bernstein 
polynomials was used. This idea was further gen- 
eralized in the works of Ref. Q where it was pro- 
posed to relate the SFs and their Mellin moments 
using the series in the Jacobi polynomials. This 
method was developed in Ref. Q and successfully 
applied in the process of the NLO fits of the DIS 
data of the BCDMS collaboration Q and of the 
old data of the CCFR collaboration of Rcl^ 
(see Ref.@). 

The Jacobi polynomial technique was already 
used at the NNLO in the process of the NS fits 
of the F2 SF data of the BCDMS collaboration 
old XF3 data of the CCFR collaboration ||l|l 
and the revised in Ref.|l5|] xFs data of the CCFR 
collaboration [|l6| (which we will denote hereafter 
as CCFR'97). 

In this work we are summarizing the definite 
results of Ref. |]l6| and are generalizing them to 
the next-to-next-to-next-to-leading order (N'^LO) 
level by means of the Fade approximants ap- 
proach. Our aim is to get the estimate of the 
value of A^IL. and as{Mz) at the N^LO and to fix 
more precisely the theoretical uncertainties of the 
values of these parameters , obtained in Ref.fl^ 
at the NNLO level. 

2. Let us define the Mellin moments for the NS 
SF xF3{x,Q^): M,f'5(g2) = x"-^F3{x,Q^)dx 
where n = 2,3,4,.... The theoretical expression 
for these moments obey the following renormal- 
ization group equation 
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where Ag — as/(47r). The renormahzation 
group functions arc defined as = PiAg) = 



X]i>o 7^s("-)^s^^ where Z^^ are the renormal- 
ization constants of the corresponding NS opera- 
tors. The solution of the renormahzation group 
equation can be presented in the following form : 



exp 



dx 



As(Ql) 



(2) 



Pix) 

where M^^{Qq) is the phenomenological quan- 
tity related to the factorization scale dependent 
factor. It can be parametrized through the parton 
distributions at fixed momentum transfer Qq as 
M,^^{Ql) = J^x^-^A(Ql)x''^Qo)(i-xYiQ'o){i + 
^{Ql)x)dx with 7 7^ 0. 

At the N'^LO the expression for the coefficient 

(n) 

function Cj^g can be presented in the follow- 
ing general form C^^(A,,) = 1 + C^^\n)As + 
C(^^(n)As -I- C^^^n)Al, while the corresponding 
expansion of the anomalous dimensions term is 



exp 



-dx 



(3) 



(A,(Q2))<^(")/2/'0x[1+P(„)^^(q2) 

+q{n)As{Q'f + r{n)A,{Q^f] 

where C'-^-'(n), C^^^(n), p{n) and q{n) are defined 
in Ref. 0. 

The coupling constant As{Q'^) can be expressed 
in terms of the inverse powers of L = \n{Q^/Aj^) 

as A^^^O ^ ^NLO + /s^^NNLO Af^o ^ 

A^NLO ^^j^fLO ^Yiere 
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Table 1. 

The results of the QCD fit of the CCFR'97 data for > SGeF^, Ql = 5GeV^ (86 
experimental points). Statistical errors are taken into account. 
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78.4 


293±29 


0.119±0.002 


65.7 


264±85 


n 1 1 7+0.006 
^■^^ ' -0.007 


N^LO 


Fade 


78.8 


307±31 


0.120±0.002 


65.7 


256±81 


'-'■^^ ' -0.007 



Notice that in our normalizations /3o, Pi, P2 
and P3 read p„ = 11 ~ 0.6667/, /?i = 102 - 
12.6667/, (32 = 1428.50 - 279.611/ + 6.01852/^ 
and (33 = 29243.0 - 6946.30/ + 405.089/^ + 
1.49931/'^ where the expression for (3^ was ob- 
tained in Ref.0. The inverse-log expansion 
for AA^ , which incorporates the information 
about the coefficient (3^ , was presented in Ref . [|l8| . 

In order to generalize the considerations of 
Refs. |p^ , |l6| to the level of the explicitly unknown 
N'^LO corrections wc will use the Fade approx- 
imations approach of Refs. |p^ , po| . This method 
provides a possibility to estimate the higher or- 
der terms of perturbation theory using the Fade 
approximants (in the case of the NLO analysis of 
the experimental data of the BCDMS and EMC 
collaborations the similar idea was proposed in 
Ref. 1^ ) . In the framework of this technique the 
values of C^'^^ (n) and r{n) could be expressed as 

Fade [1/1] : C^^\n) = [C'^^Hn)]^ / C^^\n)i7) 
r{n) = q{nf/p{n) (8) 
Pade [0/2] : C^^\n) = 2C^^\n)C^^\n) (9) 

r{n) = 2p{n)q{n) - [p{n)f (10) 
Using the following equation 

:,F^-^^{x,Q^)^^+x-{l-xf -K (11) 

£ e::-^(x)5]4")(a,/?)M,+2,.F3 iQ') , 

n=0 j=0 

one can reconstruct the SF from the correspond- 
ing Mellin moments. Here 8"''' are the Jacobi 



polynomials and a, (3 are their parameters, fixed 
by the condition of the requirement of the mini- 
mization of the error of the reconstruction of the 
SF. In our analysis we are considering the region 
6 < Nmax < 10. The target mass corrections are 
included in our fits up to order 0(Q^^)-terms. 
In the same manner as in Refs.[ p2| , p^ we are in- 
corporating in Eq.(ll) the twist-4 contribution 
h{x)/Q'^ as the additional term. 

3. The details of our NLO and NNLO fits, 
made for the case of / = 4 number of active 
fiavours, are described in Refs. [|lj,|l6|. In order 
to perform the N'^LO analysis we used the [0/2] 
Fade approximant motivated estimates of C^^^ (n) 
and r[n) terms (see Eqs.(9),(10)). It should be 
stressed that the N'^LO [1/1] Fade approximant 
description of the CCFR'97 experimental data 
turned out to be not acceptable, since it pro- 
duces rather high value of x^: ^ jnep > 2 (where 
nep = 86 is the number of the experimental points 
taken into account). The similar effect of the 
preference of the [0/2] Fade approximant anal- 
ysis over the [1/1] one was found in Ref.]^ in 
the case of the comparison of the QCD theoreti- 
cal predictions for the polarized Bjorken sum rule 
(which are closely related to the QCD predictions 
for the first moment of the xF^ SF, namely for the 
Gross-Llewellyn Smith sum rule) with the avail- 
able experimental data. 

The resuhs for A^, obtained at the LO, NLO, 

MS 

NNLO and N^LO, are presented in Table 1. Us- 
ing the LO, NLO, NNLO and N^LO variants of 
the rigorous M S-scheme matching conditions, de- 
rived in Ref. |lj] following the lines of Ref. j23| , we 
are transforming A.^^ values through the thresh- 

AI S 

old of the production of the fifth flavour M5 = mt, 
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(where rub is the 6-quark pole mass) and are ob- 
taining the related values of and as{AIz)- 

The results of the LO, NLO and NNLO ex- 
tractions of the twist-4 parameter h{x) from the 
CCFR'97 xFs data were presented in Ref.||l|]. 
Using the Pade approximations approach we also 
found the a:- form of h{x) at the N^LO. It turned 
out to be almost undistinguished from the results 
of the NNLO fits. Dew to the lack of space we are 
postponing the presentation of the N^LO results 
for h(x) for the future publication. 

In the case of the non-zero values of the twist-4 
parameter h{x) ^ the Pade motivated N'^LO 
results for A^l^ and as{Mz) turned out to be 
almost identical to the NNLO ones (provided 
the statistical error bars are taken into account, 
see Tabic 1). Moreover, the difference A^^-^'^= 



/-A (4) \N-'LO _ / A i-l) \NNLO\ 
^ MS' ^ MS' ' 

sidered as the measure of the theoretical un- 
certainties of the NNLO results, is drastically 
smaller then the NLO correction term A^-^'^= 
,(4) ^NNLO _ fA(4) ^NLOi ^he similar ten- 



which can be con- 



(Ai 



MS' ' MS' 

dency A^^-^° << A^^*^ is taking place in the 
case of the fits without twist-4 corrections. These 
observed properties indicate the reduction of the 
theoretical errors due to cutting the analyzed per- 
turbative series at the different orders. 

It is known that the inclusion of the higher- 
order perturbativc QCD corrections into the com- 
parison with the experimental data is decreasing 
the scale-scheme theoretical errors of the results 
for Ajjg and thus asiMz) (see e.g. Refs.[|4[|o)). 
Among the ways of probing the scale-scheme 
uncertainties are the scheme-invariant methods, 
namely the principle of minimal sensitivity ||25[] 
and the effective charges approach |2^], which 
is known to be identical to the scheme-invariant 
perturbation theory, developed in Refs.|27-29|. 
These methods can be used to estimate unknown 
N'^LO corrections to the definite physical quanti- 
ties [Q. Note that the estimates of Ref.H^ are 
in agreement with the results of applications of 
the Pade resummation technique (see Ref.|19|). 
Therefore, we can conclude that the application 
of the methods of the Pade approxiniants is lead- 
ing to the reduction of the scale-scheme depen- 
dence uncertainties of the values of as{Mz). We 



are presenting now the outcomes of the fits of the 
CCFR'97 experimental data for the xFs SF, ob- 
tained with considering twist-4 parameters as the 
additional free parameters of the fit: 

NLO HT free a,{Mz) = 0.125t^:j]?g(s<at) (12) 
±0m5{syst) ± 0.009{theory) 

NNLO HT free as{Mz) = 0. 117 t°Q°„°Q%stat) (13) 
±0m5{syst) ± 0m2{theory) 

where the systematic uncertainties are taken from 
the experimental analysis of Ref. |l^ and the the- 
oretical uncertainties in the results of Eq.(12) 
[eq.(13)] are estimated by the diffferences between 
the central values of the outcomes of the NNLO 
and NLO [N^LO and NNLO] fits, presented in 
Table 1, plus the arbitrariness in the application 
of the MS'-scheme matching condition of Ref. ||l^ , 
which following the considerations of Ref. ||3^ we 
estimate as Aas{Mz) = ±0.001. 

It can be seen that due to the large over- 
all number of the fitted parameters the results 
of Eqs.(12),(13) have rather large statistical un- 
certainties. It is possible to decrease their val- 
ues by fixing the concrete form of the twits- 
4 parameter h{x). For example, parametrizing 
h{x) through the infrared renormalon model of 
Ref. [Q, supported by the independent consid- 
erations of Ref. js^, we arrive to the following 
values of as{Mz), extracted in Ref.|^^ from the 
CCFR'97 data for xF^ SF: 

NLO HT of [32] asiMz) = 0.121 ± 0.002(stat) (14) 
±0m5{syst) ± Qm6{theory) 

NNLO HT of [32] as{Mz) = 0.117 ± 0.002(stat) (15) 
±0.005{syst) ± 0.003{theory) . 

Note, however, that we did not yet check the 
validity of the estimates of the theoretical error 
bars of the result of Eq.(15) with the help of the 
method of the [0/2] Pade approximant. We are 
planning to perform this analysis in another work. 
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